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A LOGICAL CALCULUS OF THE 
IDEAS IMMANENT IN NERVOUS ACTIVITY 


WARREN S. MCCULLOCH AND WALTER PITTS 


FROM THE UNIVERSITY OF ILLINOIS, COLLEGE OF MEDICINE, 
DEPARTMENT OF PSYCHIATRY AT THE ILLINOIS NEUROPSYCHIATRIC INSTITUTE, 
AND THE UNIVERSITY OF CHICAGO 


Because of the “all-or-none” character of nervous activity, neural 
events and the relations among them can be treated by means of propo- 
sitional logic. It is found that the behavior of every net can be described 
in these terms, with the addition of more complicated logical means for 
nets containing circles; and that for any logical expression satisfying 
certain conditions, one can find a net behaving in the fashion it describes. 
It is shown that many particular choices among possible neurophysiologi- 
cal assumptions are equivalent, in the sense that for every net behav- 
ing under one assumption, there exists another net which behaves un- 
der the other and gives the same results, although perhaps not in the 
same time. Various applications of the calculus are discussed. 


I. Introduction 


Theoretical neurophysiology rests on certain cardinal assump- 
tions. The nervous system is a net of neurons, each having a soma 
and an axon. Their adjunctions, or synapses, are always between the 
axon of one neuron and the soma of another. At any instant a neuron 
has some threshold, which excitation must exceed to initiate an im- 
pulse. This, except for the fact and the time of its occurrence, is de- 
termined by the neuron, not by the excitation. From the point of ex- 
citation the impulse is propagated to all parts of the neuron. The 
velocity along the axon varies directly with its diameter, from less 
than one meter per second in thin axons, which are usually short, to 
more than 150 meters per second in thick axons, which are usually 
long. The time for axonal conduction is consequently of little impor- 
tance in determining the time of arrival of impulses at points un- 
equally remote from the same source. Excitation across synapses oc- 
curs predominantly from axonal terminations to somata. It is still a 
moot point whether this depends upon irreciprocity of individual syn- 
apses or merely upon prevalent anatomical configurations. To sup- 
pose the latter requires no hypothesis ad hoc and explains known ex- 
ceptions, but any assumption as to cause is compatible with the cal- 
culus to come. No case is known in which excitation through a single 
synapse has elicited a nervous impulse in any neuron, whereas any 


115 


116 LOGICAL CALCULUS FOR NERVOUS ACTIVITY 


neuron may be excited by impulses arriving at a sufficient number of 
neighboring synapses within the period of latent addition, which lasts 
less than one quarter of a millisecond. Observed temporal summation 
of impulses at greater intervals is impossible for single neurons and 
empirically depends upon structural properties of the net. Between 
the arrival of impulses upon a neuron and its own propagated im- 
pulse there is a synaptic delay of more than half a millisecond. Dur- 
ing the first part of the nervous impulse the neuron is absolutely re- 
fractory to any stimulation. Thereafter its excitability returns rap- 
idly, in some cases reaching a value above normal from which it sinks 
again to a subnormal value, whence it returns slowly to normal. Fre- 
quent activity augments this subnormality. Such specificity as is 
possessed by nervous impulses depends solely upon their time and 
place and not on any other specificity of nervous energies. Of late 
only inhibition has been seriously adduced to contravene this thesis. 
Inhibition is the termination or prevention of the activity of one 
group of neurons by concurrent or antecedent activity of a second 
group. Until recently this could be explained on the supposition that 
previous activity of neurons of the second group might so raise the 
thresholds of internuncial neurons that they could no longer be ex- 
cited by neurons of the first group, whereas the impulses of the first 
group must sum with the impulses of these internuncials to excite the 
now inhibited neurons. Today, some inhibitions have been shown to 
consume less than one millisecond. This excludes internuncials and 
requires synapses through which impulses inhibit that neuron which 
is being stimulated by impulses through other synapses. As yet ex- 
periment has not shown whether the refractoriness is relative or ab- 
solute. We will assume the latter and demonstrate that the difference 
is immaterial to our argument. Hither variety of refractoriness can 
be accounted for in either of two ways. The “inhibitory synapse” 
may be of such a kind as to produce a substance which raises the 
threshold of the neuron, or it may be so placed that the local disturb- 
ance produced by its excitation opposes the alteration induced by the 
otherwise excitatory synapses. Inasmuch as position is already known 
to have such effects in the case of electrical stimulation, the first hy- 
pothesis is to be excluded unless and until it be substantiated, for the 
second involves no new hypothesis. We have, then, two explanations 
of inhibition based on the same general premises, differing only in 
the assumed nervous nets and, consequently, in the time required for 
inhibition. Hereafter we shall refer to such nervous nets as equiva- 
lent in the extended sense. Since we are concerned with properties 
of nets which are invariant under equivalence, we may make the 
physical assumptions which are most convenient for the calculus. 
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Many years ago one of us, by considerations impertinent to this 
argument, was led to conceive of the response of any neuron as fac- 
tually equivalent to a proposition which proposed its adequate stimu- 
lus. He therefore attempted to record the behavior of complicated 
nets in the notation of the symbolic logic of propositions. The “all- 
or-none”’ law of nervous activity is sufficient to insure that the activ- 
ity of any neuron may be represented as a proposition. Physiological 
relations existing among nervous activities correspond, of course, to 
relations among the propositions; and the utility of the representa- 
tion depends upon the identity of these relations with those of the 
logic of propositions. To each reaction of any neuron there is a corre- 
sponding assertion of a simple proposition. This, in turn, implies 
either some other simple proposition or the disjunction or the con- 
junction, with or without negation, of similar propositions, according 
to the configuration of the synapses upon and the threshold of the 
neuron in question. Two difficulties appeared. The first concerns 
facilitation and extinction, in which antecedent activity temporarily 
alters responsiveness to subsequent stimulation of one and the same 
part of the net. The second concerns learning, in which activities 
concurrent at some previous time have altered the net permanently, 
so that a stimulus which would previously have been inadequate is 
now adequate. But for nets undergoing both alterations, we can sub- 
stitute equivalent fictitious nets composed of neurons whose connec- 
tions and thresholds are unaltered. But one point must be made clear: 
neither of us conceives the formal equivalence to be a factual expla- 
nation. Per contra!—we regard facilitation and extinction as depen- 
dent upon continuous changes in threshold related to electrical and 
chemical variables, such as after-potentials and ionic concentrations ; 
and learning as an enduring change which can survive sleep, anaes- 
thesia, convlusions and coma. The importance of the formal equiva- 
lence lies in this: that the alterations actually underlying facilitation, 
extinction and learning in no way affect the conclusions which follow 
from the formal treatment of the activity of nervous nets, and the 
relations of the corresponding propositions remain those of the logic 
of propositions. 

The nervous system contains many circular paths, whose activity 
so regenerates the excitation of any participant neuron that reference 
to time past becomes indefinite, although it still implies that afferent 
activity has realized one of a certain class of configurations over time. 
Precise specification of these implications by means of recursive func- 
tions, and determination of those that can be embodied in the activity 
‘of nervous nets, completes the theory. 
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Il. The Theory: Nets Without Circles 


We shall make the following physical assumptions for our cal- 
culus. 


1. The activity of the neuron is an “all-or-none” process. 


2. <A certain fixed number of synapses must be excited within 
the period of latent addition in order to excite a neuron at any time, 
and this number is independent of previous activity and position on 
the neuron. 


9 


8. The only significant delay within the nervous system is syn- 
aptic delay. 

4. The activity of any inhibitory synapse absolutely prevents 
excitation of the neuron at that time. 


5. The structure of the net does not change with time. 


To present the theory, the most appropriate symbolism is that of 
Language II of R. Carnap (1938), augmented with various notations 
drawn from B. Russell and A. N. Whitehead (1927), including the 
Principia conventions for dots. Typographical necessity, however, 
will compel us to use the upright ‘E’ for the existential operator in- 
stead of the inverted, and an arrow (‘>’) for implication instead of 
the horseshoe. We shall also use the Carnap syntactical notations, but 
print them in boldface rather than German type; and we shall intro- 
duce a functor S, whose value for a property P is the property which 
holds of a number when P holds of its predecessor; it is defined by 
‘S(P) (t) =. P(Kx) .t = 2’)’; the brackets around its argument will 
often be omitted, in which case this is understood to be the nearest 


predicate-expression [Pr] on the right. Moreover, we shall write 
S?Pr for S(S(Pr) ), ete. 


The neurons of a given net N may be assigned designations 
‘Cy’, Ce’, +++, ‘Cn’. This done, we shall denote the property of a number, 
that a neuron ¢; fires at a time which is that number of synaptic de- 
lays from the origin of time, by ‘N’ with the numeral 7 as subscript, 
so that N;(t) asserts that ¢; fires at the time t. N; is called the action 
of c;. We shall sometimes regard the subscripted numeral of ‘N’ as 
if it belonged to the object-language, and were in a place for a func- 
toral argument, so that it might be replaced by a number-variable 
[z] and quantified; this enables us to abbreviate long but finite dis- 
junctions and conjunctions by the use of an operator. We shall em- 
ploy this locution quite generally for sequences of Pr; it may be se- 
cured formally by an obvious disjunctive definition. The predicates 
‘N,’, ‘N.’, --- , comprise the syntactical class ‘N’. 
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Let us define the peripheral afferents of N as the neurons of NX 
with no axons synapsing upon them. Let N,, --- , N, denote the ac- 
tions of such neurons and N,.;, , Np2,-:: , Nn those of the rest. Then a 
Solution of N will be a class of sentences of the form Si: Npu(2Z,) =. 
Pr;(N,, Nz ,---, Np, 2), where Pr; contains no free variable save z, 
and no descriptive symbols save the N in the argument [Arg], and 
possibly some constant sentences [sa]; and such that each S; is true 
of N. Conversely, given a Pr. (*p',, "p2,-:+, Mp, 21, 8), containing 
no free variable save those in its Arg, we shall say that it is realizable 
in the narrow sense if there exists a net N and a series of N; in it 
such that N,(z,) =. Pri(N,, No,-+:, 2,, 8a) is true of it, where sa, 
has the form N(0). We shall call it realizable in the extended sense, 
or simply realizable, if for some n S"(Pr,) (pi, --+ , Dp, 21, 8) is 
realizable in the above sense. c,; is here the realizing neuron. We 
shall say of two laws of nervous excitation which are such that every 
S which is realizable in either sense upon one supposition is also re- 
alizable, perhaps by a different net, upon the other, that they are 
equivalent assumptions, in that sense. 

The following theorems about realizability all refer to the ex- 
tended sense. In some cases, sharper theorems about narrow realiz- 
ability can be obtained; but in addition to greater complication in 
statement this were of little practical value, since our present neuro- 
physiological knowledge determines the law of excitation only to ex- 
tended equivalence, and the more precise theorems differ according 
to which possible assumption we make. Our less precise theorems, 
however, are invariant under equivalence, and are still sufficient for 
all purposes in which the exact time for impulses to pass through the 
whole net is not crucial. 

Our central problems may now be stated exactly: first, to find an 
effective method of obtaining a set of computable S constituting a 
solution of any given net; and second, to characterize the class of 
realizable S in an effective fashion. Materially stated, the problems 
are to calculate the behavior of any net, and to find a net which will 
behave in a specified way, when such a net exists. 

A net will be called cyclic if it contains a circle: i.e., if there ex- 
ists a chain ¢;, Cian, -*: Of neurons on it, each member of the chain 
synapsing upon the next, with the same beginning and end. If a set 
of its neurons ¢,, Cz ,°*:, Cp is such that its removal from WN leaves 
it without circles, and no smaller class of neurons has this property, 
the set is called a cyclic set, and its cardinality is the order of N. In 
an important sense, as we shall see, the order of a net is an index of 
the complexity of its behavior. In particular, nets of zero order have 
especially simple properties; we shall discuss them first. 
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Let us define a temporal propositional expression (a TPE), des- 
ignating a temporal propositional function (TPF), by the following 
recursion: 


1. A‘p'[z,] is a TPE, where p, is a predicate-variable. 

2. If S, and S. are TPE containing the same free individual 
variable, so are SS, , S:vS., S,.S, and Sj. © S2. 

3. Nothing else isa TPE. 


THEOREM I. 
Every net of order 0 can be solved in terms of temporal proposi- 
tional expressions. 


Let ¢; be any neuron of N with a threshold 6; > 0, and let ci, 
Cig, «+: , Cip have respectively Nir, Nie, --- , Nip excitatory synapses 
upon it. Let cj, Cj2, +++, Cjq have inhibitory synapses upon it. Let x; 
be the set of the subclasses of {Ni1 , Ni, «++ , Nip} Such that the sum of 
their members exceeds 6;. We shall then be able to write, in accord- 
ance with the assumptions mentioned above, 


Nia) =.8{ 1 © Nim(2:). SH Ni (2,)| (1) 


aéeKi Sea 


where the ‘>’ and ‘77’ are syntactical symbols for disjunctions and 
conjunctions which are finite in each case. Since an expression of 
this form can be written for each ¢c; which is not a peripheral affer- 
ent, we can, by substituting the corresponding expression in (1) for 
each N;,, or Ni; whose neuron is not a peripheral afferent, and re- 
peating the process on the result, ultimately come to an expression 
for N; in terms solely of peripherally afferent N, since N is without 
circles. Moreover, this expression will be a TPE, since obviously (1) 
is; and it follows immediately from the definition that the result of 
substituting a TPE for a constituent p(z) in a TPE is also one. 


THEOREM II. 
Every TPE ts realizable by a net of order zero. 


The functor S obviously commutes with disjunction, conjunction, 
and negation. It is obvious that the result of substituting any S;, 
realizable in the narrow sense (i.n.s.), for the p(z) in a realizable ex- 
pression S, is itself realizable i.n.s.; one constructs the realizing net by 
replacing the peripheral afferents in the net for S, by the realizing neu- 
rons in the nets for the S;. The one neuron net realizes p, (23) tds. 
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and Figure 1-a shows a net that realizes Sp,(z,) and hence SS, , i.n.s., 
if S, can be realized in.s. Now if S. and S3 are realizable then 
S"S, and S"S; are realizable i.n.s., for suitable m and n. Hence so 
are S”*"S, and S™"S; . Now the nets of Figures 1b, c and d respectively 
realize S (p1 (21) V Po(21)), S (ps (21). Do (21) ), and S(p, (21). © Po (2:) ) 
i.n.s. Hence S™*™ (S, v S.), S™™? (S, .S.), and S™ (S,.© S,) are 
realizable i.n.s. Therefore S, v S. S, . S. S: . © Se are realizable if 
S, and S, are. By complete induction, all TPE are realizable. In this 
way all nets may be regarded as built out of the fundamental elements 
of Figures la, b, c, d, precisely as the temporal propositional expres- 
sions are generated out of the operations of precession, disjunction, 
conjunction, and conjoined negation. In particular, corresponding to 
any description of state, or distribution of the values true and false 
for the actions of all the neurons of a net save that which makes 
them all false, a single neuron is constructible whose firing is a neces- 
sary and sufficient condition for the validity of that description. More- 
over, there is always an indefinite number of topologically different 
nets realizing any TPE. 


THEOREM III. 


Let there be given a complex sentence S, built up in any manner out 
of elementary sentences of the form p(z, — zz) where 2Z is any nu- 
meral, by any of the propositional connections: negation, disjunction, 
conjunction, implication, and equivalence. Then S, is a TPE and 
only if it is false when its constituent p(z, — zz) are all assumed 
false—i.e., replaced by false sentences — or that the last line in tts 
truth-table contains an ‘F’,—or there is no term in its Hilbert dis- 
junctive normal form composed exclusively of negated terms. 


These latter three conditions are of course equivalent (Hilbert 
and Ackermann, 1938). We see by induction that the first of them is 
necessary, since p(z, — zz) becomes false when it is replaced by a 
false sentence, and S, v S2, Si . S, and S, . ~ S, are all false if 
both their constituents are. We see that the last condition is sufficient 
by remarking that a disjunction is a TPE when its constituents are, 
and that any term 


SioSs se. Sait ORT is ern 
can be written as 
(SS: 8. Sa) FSV Sime Vie VSS), 


which is clearly a TPE. 


The method of the last theorems does in fact provide a very con- 
venient and workable procedure for constructing nervous nets to or- 
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der, for those cases where there is no reference to events indefinitely 
far in the past in the specification of the conditions. By way of ex- 
ample, we may consider the case of heat produced by a transient 
cooling. 

If a cold object is held to the skin for a moment and removed, 
a sensation of heat will be felt; if it is applied for a longer time, the 
sensation will be only of cold, with no preliminary warmth, however 
transient. It is known that one cutaneous receptor is affected by heat, 
and another by cold. If we let N, and Nz be the actions of the respec- 
tive receptors and N; and N, of neurons whose activity implies a sen- 
sation of heat and cold, our requirements may be written as 


N,(t).2-= =N; (1) Viet 3) et) 
N,(t).. = .Ne(t—2) . Ne (t—1) 


lll Ih 


where we suppose for simplicity that the required persistence in the 
sensation of cold is say two synaptic delays, compared with one for 
that of heat. These conditions clearly fall under Theorem III. A net 
may consequently be constructed to realize them, by the method of 
Theorem II. We begin by writing them in a fashion which exhibits 
them as built out of their constituents by the operations realized in 
Figures la, b, c, d: i.e., in the form 


N,(t) .=.S{N,(t) v SL(SNa(t)) «© No(t)]}- 
N,(t) .=.S{[SNo(t)] . No(t)}. 


First we construct a net for the function enclosed in the greatest 
number of brackets and proceed outward; in this case we run a net of 
the form shown in Figure la from c, to some neuron ¢,, say, so that 


Na(t) .=.SN2(¢). 


Next introduce two nets of the forms 1c and 1d, both running from 
C, and C, , and ending respectively at c, and say c,. Then 


N,(t) .=.S[Na(t) .No(t)].=.S[(SN2(t)) . No(t)]. 
N,z(t) .=.S[Na(t) .~ No(t)] =. S[(SN2(t)) . © No(t)]. 


Finally, run a net of the form 1b from c, and ¢, to ¢; , and derive 


N;(t) .=.S[Ni(t) vN,(t)] 
.=.S{N,(t) vS[(SN2(t)] . © No(t)}. 
These expressions for N,(t) and N,(t) are the ones desired ; and the 
realizing net in toto is shown in Figure le. 
This illusion makes very clear the dependence of the correspon- 
dence between perception and the “external world” upon the specific 
structural properties of the intervening nervous net. The same illu- 


Ill 
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sion, of course, could also have been produced under various other 
assumptions about the behavior of the cutaneous receptors, with cor- 
respondingly different nets. 

We shall now consider some theorems of equivalence: i.e., theo- 
rems which demonstrate the essential identity, save for time, of vari- 
ous alternative laws of nervous excitation. Let us first discuss the 
case of relative inhibition. By this we mean the supposition that the 
firing of an inhibitory synapse does not absolutely prevent the firing 
of the neuron, but merely raises its threshold, so that a greater num- 
ber of excitatory synapses must fire concurrently to fire it than would 
otherwise be needed. We may suppose, losing no generality, that the 
increase in threshold is unity for the firing of each such synapse; we 
then have the theorem: 


THEOREM IV. 


Relative and absolute inhibition are equivalent in the extended 
sense. 


We may write out a law of nervous excitation after the fashion 
of (1), but employing the assumption of relative inhibition instead; 
inspection then shows that this expression is a TPE. An example of 
the replacement of relative inhibition by absolute is given by Figure 
1f. The reverse replacement is even easier; we give the inhibitory 
axons afferent to c; any sufficiently large number of inhibitory syn- 
apses apiece. 

Second, we consider the case of extinction. We may write this 
in the form of a variation in the threshold 6;; after the neuron c; has 
fired; to the nearest integer—and only to this approximation is the 
variation in threshold significant in natural forms of excitation—this 
may be written as a sequence 6; + b; for 7 synaptic delays after firing, 
where b; = 0 for 7 large enough, say 7 = M or greater. We may then 
state 


THEOREM V. 
Extinction is equivalent to absolute inhibition. 


For, assuming relative inhibition to hold for the moment, we 
need merely run M circuits T;, T,, --- Ty containing respectively 1, 
2,---, M neurons, such that the firing of each link in any is sufficient 
to fire the next, from the neuron c; back to it, where the end of the 
circuit T; has just b; inhibitory synapses upon ¢;. It is evident that 
this will produce the desired results. The reverse substitution may be 
accomplished by the diagram of Figure 1g. From the transitivity of 
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replacement, we infer the theorem. To this group of theorems also 
belongs the well-known 


THEOREM VI. 


Facilitation and temporal summation may be replaced by spatial 
summation. 


This is obvious: one need merely introduce a suitable sequence 
of delaying chains, of increasing numbers of synapses, between the 
exciting cell and the neuron whereon temporal summation is desired 
to hold. The assumption of spatial summation will then give the re- 
quired results. See e.g. Figure 1h. This procedure had application in 
showing that the observed temporal summation in gross nets does not 
imply such a mechanism in the interaction of individual neurons. 

The phenomena of learning, which are of a character persist- 
ing over most physiological changes in nervous activity, seem to re- 
quire the possibility of permanent alterations in the structure of nets. 
The simplest such alteration is the formation of new synapses or 
equivalent local depressions of threshold. We suppose that some ax- 
onal terminations cannot at first excite the succeeding neuron; but if 
at any time the neuron fires, and the axonal terminations are simul- 
taneously excited, they become synapses of the ordinary kind, hence- 
forth capable of exciting the neuron. The loss of an inhibitory syn- 
apse gives an entirely equivalent result. We shall then have 


THEOREM VII. 
Alterable synapses can be replaced by circles. 


This is accomplished by the method of Figure li. It is also to be 
remarked that a neuron which becomes and remains spontaneously 
active can likewise be replaced by a circle, which is set into activity 
by a peripheral afferent when the activity commences, and inhibited 
by one when it ceases. 


III. The Theory: Nets with Circles. 


The treatment of nets which do not satisfy our previous assump- 
tion of freedom from circles is very much more difficult than that case. 
This is largely a consequence of the possibility that activity may be 
set up in a circuit and continue reverberating around it for an in- 
definite period of time, so that the realizable Pr may involve refer- 
ence to past events of an indefinite degree of remoteness. Consider 
such a net N, say of order p, and let c,, c2,+++ , ¢» be a cyclic set of 
neurons of XN. It is first of all clear from the definition that every N, 
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of N can be expressed as a TPE, of N,, Nz,---, N, and the absolute 
afferents; the solution of N involves then only the determination of 
expressions for the cyclic set. This done, we shall derive a set of ex- 
pressions [A]: 


N; (21) «=. Pri[S Ny (21), 8S" No (21), +--+, S"* Np(z:)], (2) 


where Pr; also involves peripheral afferents. Now if m is the least 
common multiple of the n;;, we shall, by substituting their equiva- 
lents according to (2) in (3) for the N;, and repeating this process 
often enough on the result, obtain S of the form 


Ni (24) .=. Pri[S” Ni (21), S" Nz (21), ++, 5" N,(2)]. (3) 


These expressions may be written in the Hilbert disjunctive normal 
form as 


Ni(4) .=.> Sall 8S" N;(z,) I © 8S" N;(z:), for suitablex, ( 4) 
ete jeK JEBa 

where S, is a TPE of the absolute afferents of N. There exist some 
2° different sentences formed out of the p N; by conjoining to the con- 
junction of some set of them the conjunction of the negations of the 
rest. Denumerating these by X1(21) , X2(z1),--- , X,, (2) , we may, 
by use of the expressions (4), arrive at an equipollent set of equations 
of the form 


Xi(z). =. 3 Pry (z,). 9X2). (5) 


j= 


Now we import the subscripted numerals 7,7 into the object-language: 
i.e., define Pr, and Pr. such that Pr,(22,,21) . =. Xi (2) and Pr2(22, 
229, 2) . =. Pri;(z,) are provable whenever zz, and zz, denote 7 and 


4 respectively. 
Then we may rewrite (5) as 
(21) ZZ) : Pry (21 ; Zs) 
ies) 2a Pra (2) ze), 2a ea) ls Co eee) (6) 
where zz, denotes n and zz, denotes 2”. By repeated substitution we 
arrive at an expression 
(2,) 22 1 Pri (21, 22Zn 222) «=. (E22) 22 (E23) 22)... (Zn) 2% - 
Pr (2: 5 2; Zin @22 — 1)). Pra (2a, 22 , 22n (2% —1))....- (7) 
Pr. (Zn) Zn 50). Pr(2Zn, 0), for any numeral zz, which denotes s . 
This is easily shown by induction to be equipollent to 
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(2) 22, SEPIA 220220) =: (EF) (22) ZZ. — 1 f (22 22n) 
< 22). f (Z2n 220) = 2%. Pr, (f (2Zn (22 + 1)), (8) 
f (z2n Z2)) . Pr (Ff (0) , 0) 


and since this is the case for all zz. , it is also true that 
(22) (4) 225 PTs (25 22) oe (EF) (22) (2, — 1) . fF (22) 


Songz, f (2{) = 2, L228. (Pref (aa DF 2s) Pee. (9) 
Pr,[f (ves (2, , 2Zn) ), reS (Zs, ZZn) J, 


where zz, denotes n , res(7,S) is the residue of 7 mod s and zz, denotes 
2°. This may be written in a less exact way as 


Ni(t) .=. (E$) («)t —1.6(x) $2. g(t) =. 
P[o(% +1), $(%) .Now)(0)], | 


where x and ¢ are also assumed divisible by n, and Pr. denotes P. 
From the preceding remarks we shall have 


THEOREM VIII. 


The expression (9) for neurons of the cyclic set of a net N together 
with certain TPE expressing the actions of other neurons in terms 
of them, constitute a solution of N. 


Consider now the question of the realizability of a set of S;. A 
‘first necessary condition, demonstrable by an easy induction, is that 


(22) 21 - Pr (22) = Po(Z2) . >. Si = Si | - | (10) 


should be true, with similar statements for the other free p in Si: i.e., 
no nervous net can take account of future peripheral afferents. Any 


S; satisfying this requirement can be replaced by an equipollent S of 
the form 


(EF) (22) 2, (23) 22, * Fe Praj 
2 F (21, 22,23) =1.=. p, (z) (11) 
where zz, denotes p , by defining 
A 
PP mi = FL (2) (22) 21 (23) 229-2 - F (21, Za.9 Sg) = 0.V.F (ey, zo Ze) 
cao Ney a -eree tee Oe bts) Peat — 8 P,, (22)2 >= Ss); 
Consider now these series of classes a; , for which 
Ni(t) : =: (E¢) (x)t(m)q: geai: Nn (x) .=. O(b 2m) =I 
[t=q+1,---,M] O23 
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holds for some net. These will be called prehensible classes. Let us 
define the Boolean ring generated by a class of classes x as the aggre- 
gate of the classes which can be formed from members of « by re- 
peated application of the logical operations; i.e., we put 


R(x) == pAl (a, 8) 20 ex 
aera, psi. —a,e. Pav Bes]% 


We shall also define 


R(x) =. R(x) — pi — “en, 
Reta il eep awake > aed o>. a a Re eS Rae] 
Re(x) =Re(x) — 1p — “en, 
and 
o(y,t) =$[(m) .g(t+1,t,m) =y(m)]. 


The class R.(x) is formed from « in analogy with R(x), but by re- 
peated application not only of the logical operations but also of that 
which replaces a class of properties P «a by S(P)e S “a. We shall 
then have the 


LEMMA 


Pr,(D; 5 P25°**»Pm,2:) isa TPE if and only if 


(Zz) (Pi Bers Pm) (EDins) * Dm € Re ({D: » Doses Pm}) (13) 
Diver (21) = Pr,(Pi, D2, ***» Pm» 21) 
is true; and it is a TPE not involving ‘S’ if and only if this holds when 


‘R,’ is replaced by ‘R’, and we then obtain 


THEOREM IX. 


A series of classes a, , dz, ++: ds 1s a Series of prehensible classes tf 
and only if 


(Em) (En) (p)n(i) (py) +. (2) myp(x) =Ovy(z) =1: >: (EB) 
(Ey)m .y(y) =0.Be Rly (Ht) .y=ai)) .v. (am. 
(2) =0.Be Rly ((Et) .y=ai)] : (t)($) : pea. 
a(o,nt+p).—>. (Ef) .feB. (w)m(x)t—-1. 
p(n(t + 1) +p,nx+p,w)=—f(nt+p,nx+p,w). 


(14) 
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The proof here follows directly from the lemma. The condition is 
necessary, since every net for which an expression of the form (4) 
can be written obviously verifies it, the y’s being the characteristic 
functions of the S. and the # for each w being the class whose desig- 
nation has the form 7 Pr; IJ Pr;, where Pr, denotes o, for all #. Con- 

iea g&Ba 
versely, we may write an expression of the form (4) for a net N 
fulfilling prehensible classes satisfying (14) by putting for the Pra 
Pr denoting the y’s, and a Pr, written in the analogue for classes of 
the disjunctive normal form, and denoting the a corresponding to that 
wy, conjoined to it. Since every S of the form (4) is clearly realizable, 
we have the theorem. 

It is of some interest to consider the extent to which we can by 
knowledge of the present determine the whole past of various special 
nets: i.e., when we may construct a net the firing of the cyclic set of 
whose neurons requires the peripheral afferents to have had a set of 
past values specified by given functions ¢;. In this case the classes 
a; of the last theorem reduced to unit classes; and the condition may 
be transformed into 


(Em,n) (p)n(t, yp) (E7) :. (a) m: yp(“%) =0 Nv. yp(x)=1: 
diea(y,nt+p):>: (w)m(a2)t—-1.¢4:(n(t+1) 
+p, natp,w)=d(nt+p,ne+p,w) :. 

(u,v) (w)m.di(n(ut+1)+p,nut+p,w) 
=di(n(v+1)+p,nv+o0,w). 


On account of limitations of space, we have presented the above 
argument very sketchily; we propose to expand it and certain of its 
implications in a further publication. 

The condition of the last theorem is fairly simple in principle, 
though not in detail; its application to practical cases would, however, 
require the exploration of some 2” classes of functions, namely the 
members of R({a, ,--- , as}). Since each of these is a possible B of 
Theorem IX, this result cannot be sharpened. But we may obtain a 
sufficient condition for the realizability of an S§ which is very easily 
applicable and probably covers most practical purposes. This is given 
by 


THEOREM X. 


Let us define a set K of S by the following recursion: 


1. Any TPE and any TPE whose arguments have been replaced 
by members of K belong to K; 


2. IfPr,(z,) is a member of K , then (22)2, . Pri(z.), (E22) 2 . 
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Pr,(22), and Cmn(z1) . 8 belong to it, where C,,, denotes the property 
of being congruent to m modulo n,m <n. 


3. The set K has no further members 
Then every member of K is realizable. 
For, if Pr,(z,) is realizable, nervous nets for which 


Ni (21) .=. Pry (21) «SN; (2:1) 
Nj (2,) .=. Pr, (21) v SN; (2:1) 


are the expressions of equation (4), realize (z.)z, . Pr,(z.) and 
(HE z2)z, . Pri(z.) respectively; and a simple circuit, c,, Co ,--- , Cr» 
of x links, each sufficient to excite the next, gives an expression 


Nn (21) -=.Ni(0) «Cin 


for the last form. By induction we derive the theorem. 

One more thing is to be remarked in conclusion. It is easily 
shown: first, that every net, if furnished with a tape, scanners con- 
nected to afferents, and suitable efferents to perform the necessary 
motor-operations, can compute only such numbers as can a Turing 
machine; second, that each of the latter numbers can be computed by 
such a net; and that nets with circles can be computed by such a net; 
and that nets with circles can compute, without scanners and a tape, 
some of the numbers the machine can, but no others, and not all of 
them. This is of interest as affording a psychological justification of 
the Turing definition of computability and its equivalents, Church’s 
A — definability and Kleene’s primitive recursiveness: If any number 
can be computed by an organism, it is computable by these definitions, 
and conversely. 


IV. Consequences 


Causality, which requires description of states and a law of nec- 
essary connection relating them, has appeared in several forms in 
several sciences, but never, except in statistics, has it been as irre- 
ciprocal as in this theory. Specification for any one time of afferent 
stimulation and of the activity of all constituent neurons, each an 
“ajl-or-none” affair, determines the state. Specification of the ner- 
vous net provides the law of necessary connection whereby one can 
compute from the description of any state that of the succeeding 
state, but the inclusion of disjunctive relations prevents complete de- 
termination of the one before. Moreover, the regenerative activity 
of constituent circles renders reference indefinite as to time past. 
Thus our knowledge of the world, including ourselves, is incomplete 


AAA AE 
fs Ml Ae 
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as to space and indefinite as to time. This ignorance, implicit in all 
our brains, is the counterpart of the abstraction which renders our 
knowledge useful. The role of brains in determining the epistemic 
relations of our theories to our observations and of these to the facts 
is all too clear, for it is apparent that every idea and every sensation 
is realized by activity within that net, and by no such activity are the 
actual afferents fully determined. 

There is no theory we may hold and no observation we can make 
that will retain so much as its old defective reference to the facts if 
the net be altered. Tinitus, paraesthesias, hallucinations, delusions, 
confusions and disorientations intervene. Thus empiry confirms that 
if our nets are undefined, our facts are undefined, and to the “real’’ 
we can attribute not so much as one quality or “form.” With determi- 
nation of the net, the unknowable object of knowledge, the “thing in 
itself,” ceases to be unknowable. 

To psychology, however defined, specification of the net would 
contribute all that could be achieved in that field—even if the analysis 
were pushed to ultimate psychic units or ‘“‘psychons,” for a psychon 
can be no less than the activity of a single neuron. Since that activity 
is inherently propositional, all psychic events have an intentional, or 
“semiotic,” character. The “all-or-none” law of these activities, and 
the conformity of their relations to those of the logic of propositions, 
insure that the relations of psychons are those of the two-valued logic 
of propositions. Thus in psychology, introspective, behavioristic or 
physiological, the fundamental relations are those of two-valued logic. 


EXPRESSION FOR THE FIGURES 


In the figure the neuron c; is always marked with the numeral 7 upon the 
body of the cell, and the corresponding action is denoted by ‘N’ with 7 as sub- 
script, as in the text. 


Figurela N,(t) . 
Figure lb N,(¢) . 
Figurelec N,(t) . 
Figureld WN,(¢) . 
Figure le N,(€): 

N,(¢) . 
Figurelf N,(¢): 


Wee 1 

NOG 1) WN, (t— 1) 

N= 1) NA (b= 1) 

NOG © Nh) 

: N,(t —1) .v.N,(t— 8). © N,(t — 2) 

IN, (@— 2). Nz (t— 1) 

: ON, (t — 1). No(t — 1) vNg(t — 1). v- Ny(t — 1). 
N,(¢ — 1). N,(f— 1) 

: © N,(t — 2). Na (t — 2) vy N,(t — 2) .v. N,(t — 2). 
N,(t — 2). NN, — 2 

Nj (t — 2). © N,(t — 8) 

-N,(t—1).N,(t — 2) 

:N,(t —1).v-N,(t —1).(Ex)t —1.N, (x) .N, (2) 


eM aU 


N,(t) : 


Figurelg WN,(¢) . 
Figure ih N,(¢) . 
Figureli N,(t): 
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Hence arise constructional solutions of holistic problems involving 
the differentiated continuum of sense awareness and the normative, 
perfective and resolvent properties of perception and execution. From 
the irreciprocity of causality it follows that even if the net be known, 
though we may predict future from present activities, we can deduce 
neither afferent from central, nor central from efferent, nor past from 
present activities—conclusions which are reinforced by the contra- 
dictory testimony of eye-witnesses, by the difficulty of diagnosing 
differentially the organically diseased, the hysteric and the malingerer, 
and by comparing one’s own memories or recollections with his con- 
temporaneous records. Moreover, systems which so respond to the 
difference between afferents to a regenerative net and certain activity 
within that net, as to reduce the difference, exhibit purposive beha- 
vior; and organisms are known to possess many such systems, sub- 
serving homeostasis, appetition and attention. Thus both the formal 
and the final aspects of that activity which we are wont to call mental 
are rigorously deduceable from present neurophysiology. The psy- 
chiatrist may take comfort. from the obvious conclusion concerning 
causality—that, for prognosis, history is never necessary. He can 
take little from the equally valid conclusion that his observables are 
explicable only in terms of nervous activities which, until recently, 
have been beyond his ken. The crux of this ignorance is that infer- 
ence from any sample of overt behavior to nervous nets is not unique, 
whereas, of imaginable nets, only one in fact exists, and may, at any 
moment, exhibit some unpredictable activity. Certainly for the psy- 
chiatrist it is more to the point that in such systems “Mind” no longer 
“goes more ghostly than a ghost.” Instead, diseased mentality can 
be understood without loss of scope or rigor, in the scientific terms of 
neurophysiology. For neurology, the theory sharpens the distinction 
between nets necessary or merely sufficient for given activities, and 
so clarifies the relations of disturbed structure to disturbed function. 
In its own domain the difference between equivalent nets and nets 
equivalent in the narrow sense indicates the appropriate use and im- 
portance of temporal studies of nervous activity: and to mathemati- 
cal biophysics the theory contributes a tool for rigorous symbolic 
treatment of known nets and an easy method of constructing hypo- 
thetical nets of required properties. 
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A formal method is derived for converting logical relations among 
the actions of neurons in a net into statistical relations among the fre- 
quencies of their impulses. 


Consider the neuron c, upon which ¢ , -:: , Cp have excitatory, 
and Cpr , *** » Cpiq, inhibitory synapses. Let 6 be the period of latent 
addition, so that c, is excited at the time ¢ if and only if the number 
of impulses along c, ,--- , ¢) which have concurred within an interval 
of duration 6 about ¢t — 1 exceeds 6, , and none has occurred within 6 
along any of Cy: , ++: , Cpq. Suppose that the sequences of impulses 
along ¢, , --- , ¢, are Statistically independent, and let »; be the mean 
frequency of impulses along c;. Then the mean proportion of time 
that exactly r excitatory impulses arrive upon c, within an interval 
of duration 6 is given by 


tr=p to=p 9 ti=p j=r ’ 
Pai Macciojay sae py EAN ae res (1) 
tp=ty—41 49-7341) ~43=1 J=1 
Consequently, the mean proportion of time A that the number of im- 
pulses concurring upon ¢, within 6 exceeds 6, is obtained by summing 
(1) from r = 6, tor =. The intervals of duration 6. within which 
impulses from any of Cpi1, -++ , Cpq OCCUY are intervals when c, is inex- 
citable. The mean proportion of time when this is not the case is 
given by 
k=p+q 


Tiel — os). (2) 


k=p+1 


The frequency v, is then the product of (2) by A, or 


k=p+q T=p ir=p t=p t3=p j=r 
ors Ol L a0 Mie, coh ei seh IT 6". (3) 
k=p+1 1=0n jpmtpatl ge t41 44=1 f=1 


If we compare equation (3) with expression (1) of the preceding 
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paper (McCulloch and Pitts, 1943), whose notation we shall use, and 
apply a straightforward inductive argument, we obtain the following 


THEOREM. 
Let N be a net of order zero, or, more generally, one for which 
N;(z) =Sit=1 an , S| 


is a solution wherein the S; fulfill the conditions of Theorem X of the 
preceding paper. Let v; be the mean frequency of impulses in c; , and 
let the expression A, be generated out of S; by the following rules: 


(1) Replace each N of S, by ‘5 v’ with the same subscript. 
(2) Replace every v by ‘ +’, every - by ‘ X’ and every © by 


4 ae 
Z2=Z1 


(3) Replace the operators (z,) and (Ez;) respectively by I 
zi1=0 


and > wherever they occur. 
z3=0 
(4) Replace every occurrence of a predicate Cnn by a symbol for 
the function fmn(t) which is defined for all natural numbers t as unity 
when t = m(mod n), and otherwise zero. 
Then the frequencies of impulses in the ¢; are given in terms of 
those of the peripheral afferents by the equations 


6¥,—A,[i=—1,---,8]. (4) 


The correspondence expressed by this theorem is exactly that of 
Boole between the algebra of logic and that of probability. It con- 
nects the logical calculus of the preceding paper with previous treat- 
ments of the activity of nervous nets in mathematical biophysics and 
with quantitatively measurable psychological phenomena. For these 
phenomena we can construct hypothetical nets by the powerful meth- 
ods of the preceding calculus. The theorem then enables us to deter- 
mine specific predictions from the quantitative characters of the 
stimulus to those of the response. These predictions can be compared 
with observations and, if necessary, the nets be altered until the con- 
sequent predictions are verified. 

But this procedure leads to error whenever the assumptions lead- 
ing to equation (8) are not fulfilled. This will be the case if the fre- 
quencies are too great, but the limit is many times the maximum ob- 
served. When the frequencies are too small, the statistical treatment, 
though valid, is of little help, but we may here conveniently use the 
logical calculus directly. The same obtains in microscopic physiologi- 
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cal analysis. The most important exceptions are those which arise 
from physiological factors that synchronize activities and thus re- 
strict the domain of the validity of the assumption of statistical in- 
dependences required in the derivation of equation (8). 
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The theory introduces two variables ¢ and ¥. The first represents 
the intensity of emotion, the second measures the intensity of activity. 
A set of integrodifferential equations is assumed to govern the varia- 
tion of ¢ and ¥ with respect to time. Since for increasing values of ¢ the 
conduct of the organism varies from great impassivity through a normal 
level of feeling to extremes of a circular depression or catatonic excite- 
ment; whereas an increase of y results in a transition from stupor to 
manic excitement, the solutions of the equations represent quantitative 
specifications of different psychotic states. 


In the following discussion we shall develop a mathematical the- 
ory of the group of mental disorders which may be characterized in 
the following fashion: That the characteristic course of the disease 
may be essentially described in terms of the vicissitudes of two vari- 
ables, the first representing the level of feeling, affect, or emotion in 
the organism, and the other the level of activity or conation. This 
group we consider to comprise the circular insanities, the reactive 
psychoses, and the catatonia of Kahlbaum. It may also be supposed 
to include affective disorders superimposed upon psychoses of an- 
other type, and perhaps also, with a more specific interpretation of 
the determining variables, certain forms of neurosis; but we shall 
not consider these latter cases in detail. 

We shall denote the fundamental quantities of our theory, re- 
garded as functions of the time t , by the symbols ¢(t) and p(t). For 
increasing values of ¢(¢) the conduct of the organism will vary from 
great impassivity to a normal level of feeling, then to normal strong 
emotion, and finally to the extremes of a circular depression or cata- 
tonic excitement. (We do not here take account of the quality of the 
emotion). If w(t) rises indefinitely, the organism will pass from the 
greatest stupor and inactivity to normal conation, and ultimately to 
manic excitement or an idiomotor seizure. We shall measure these 
variables from an origin placed at the resting values normal for the 
organism, so that positive values of ¢(t) and w(t) correspond to 
supernormal affective and conative levels, and negative values to sub- 
normal! levels. 
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Consider the change A ¢ of $(t) during a small interval of time 
t,t + At. This quantity may be supposed to arise as the sum of 
three contributions, say 4: ¢, 42 ¢, As ¢ The first represents the 
influence of the contemporary environment at the time in question ; 
for us this can be taken only as given, and specified as an empirical 
function M(t). The second, A. ¢, represents the effect of the previ- 
ous experience of the organism; we may arrive at an expression for 
it in the following fashion. 

The quantity 4, ¢(t) for any subject is an average over a large 
number of terms, each of which represents the affectivity of the re- 
action he habitually makes to some particular kind of situation. Let 
us denote one of these terms by Az ¢a(t), and consider the contribu- 
tion to it of the events occurring during some very short interval of 
time in the past, say that from 7 until 7 + dy, where0 =7<t. De- 
note this contribution by Ay Az ¢a(t) « 

We observe first that 4, 42 ¢4(t) depends upon the degree to 
which the organism is concerned with situations of the given type at 
the time 7; times when he is little or not at all concerned with them 
will have a negligible direct effect upon his future reactions to them, 
whereas periods of very great concern will usually have a very large 
effect. Denoting this concern by a function Q.:(7), we shall embody 
it as a factor in the expression for A, As ¢4, as is simplest to satisfy 
these conditions. Suppose then that Q(7) > 0, and let the subject’s 
concern at time 7 with the given kind of situation have been attended 
with emotion, so that ¢(7) is large. By the principles of learning, 
this occasion will cause recurrences of this type of matter to be treated 
with greater emotion than otherwise, or to manifest a higher value 
of ¢; and this will be the larger, as $(7) is greater. But if the sub- 
ject’s reactions at the time 7 have been of primarily conative type, 
with a high value for w(y7) and a correspondingly reduced emotion, 
the recurrences of the same kind of situation will also exhibit a small- 
er ¢. This is to say, that A, A. ¢4(t) increases with the value of 6(n) 
and decreases with y(7), while the analogously defined Ay As walt) 
behaves in an inverse fashion. 

We shall suppose that this relation is a simple proportionality, or 


An Az oa(t) = B’ Qa(n) [6 (n) — es: p(n) ] 
Ay As pat) =f’ Quln) [p(q) — 206 (y)] (1) 


in which f’, «,, and «, are suitable constants of proportionality. 

But the strength of any contribution to A, ¢a(t) from the past 
will in general be less in an amount dependent upon its remoteness in 
time, by ordinary forgetting; so that we must multiply the left of 
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equation (1) by some decreasing function of the distance t — n into 
the past. For this we shall select the function e-”, which seems to be 
a fair representation of the “law of forgetting”, according to Ebbing- 
haus and others; the undetermined coefficient u here allows us to re- 
serve to psychoanalysts and other people who do not consider that an 
affecting experience ever becomes less efficient, the option of setting 
u = 0. We shall see, however, that this assumption requires greatly 
increased strength in the homeostatic mechanism of the organism to 
avoid continuously disturbed behavior. 

The total magnitude of A, ¢(t) will be obtained by adding to- 
gether all the quantities e#'” A, A, d4(t) corresponding to a set of 
past intervals of time which covers (0, t) completely, and then sum- 
ming quantities of this kind for each type of situation: Being linear 
they superpose, and if we allow the intervals into which we have di- 
vided a past to increase without limit in number, and their lengths 
to become uniformly indefinitely small, the error we have made in 
supposing that ¢ and y» do not change in the interior of each interval 
approaches zero, and the sum will be transformed into an integral. 
We derive 


A, p(t) =B'S Q(n)er"? [o(y) — er v(y) I] ay, (2) 


together with a corresponding expression for A, y(t), where Q(t) is 
the same kind of average of the Q,’s that $(t) is of the ¢,’s. This 
Q(t) is in principle determinable empirically, in terms of the intui- 
tive characterization we have given it above, but not easily so in 
fact. We may approximate it in somewhat the following fashion. On 
the average, those matters will tend to be of concern to a subject 
which excite his feelings, emotions, activity; and they will be of the 
more concern to him, the more they do so. This consideration will 
perhaps make it plausible to put Q = y ¢ into equation (2) and its 
analogue for y , where y is a constant of proportionality. 

The last part of A ¢, the term A; $(t), arises out of general 
homeostatic mechanisms of the organism. These tend to keep the val- 
ues of ¢ and y close to the normal levels: whenever a deviation occurs 
from these levels, it brings into play restoring forces whose strength 
increases with the extent of the deviation and tends to reduce it. We 
may represent this component of A ¢ and A y by a function f , so that 
As 6 = — f(¢), 4s y = — f(y), where f is properly monotone and 
never negative. If we now pass to the limit as 4 ¢ becomes small, and 
set 6 = f’, we shall obtain the equations 
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TF pf! sinew ta (a) —e, p(n] dyn+M(t) — flo], 


(3) 
= py b (yg) eum [p(n) — &26(n)] dn +M(t) — f[o(t)]. 
These equations form the basis of our theory, and, with the addition 
of a few subsidiary hypotheses about organic interventions, every- 
thing is to be deduced from them. In the present study we shall in 
general treat M and M , which represent the external fortunes or mis- 
fortunes of the organism, as if they consisted of a series of impulsive 
shocks between which they are relatively negligible. The “small pains 
and troubles of daily existence,” in Schopenhauer’s phrase, which 
most reasonably go into this term, may be taken account of in the de- 
termination of the normal values whence we measure ¢ and y. 

Our first problem is to determine the simplest form for the re- 
storing function f(x) which will prevent the functions ¢ and wy from 
increasing or decreasing without bound and assuming values which 
are physiologically meaningless. Considering f(z) to be a polynomial, 
we see that its leading term must be an odd power of x , in order that 
the homeostasis may work against both positively and negatively ab- 
normal values. If it is in addition not of the first degree, we may 
demonstrate its adequacy to limit ¢ and y in the following way. 

Multiply the equations (3) by e#' throughout, differentiate, and 
cancel the exponential factor. We have then 


¢°=—[uet+f(e)] ¢ +h e—fpeadoy— uf(¢) 
v= Te P@)I¥ F8eb= Fee Rae 


with the initial conditions ¢’(0) = — f[¢(0)], w'(0) = — f[w(0)], 
where we have dropped the terms in M and M, since the variables 
cannot become infinite during the integrable impulses of which M 
and M consist, and at other times they obey the equations (4). Now 
time may be divided into two sorts of interval: those where |¢! = |y|, 
and those where the reverse inequality holds. Considering intervals 
of the former type, we find, upon rearranging the first of equations 


(4), 


i sake Sas a $” uf (¢) 
ut f'(¢) uetf(¢) wtf (¢) 
$" (5) 
='0(1) — O(¢) — 


ut f(g)’ 
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as  — o through intervals of the given type, under the given assump- 
tions. Now we may distinguish these kinds of intervals of time. As 
¢ increases through intervals of the first kind, we shall have ¢” = 
O[f' (4) ]; in this case, by (6) ¢' is monotonically decreasing without 
bound for sufficiently large ¢, and hence will ultimately become nega- 
tive. As ¢ increases through intervals of the second kind, $" is nega- 
tive and of a greater order than f’($): ¢” = — 4 f' (¢), for some 
4 >0. Since f’(¢) increases without bound with its argument, in in- 
tervals of this type ¢’ will also ultimately decrease monotonically and 
without bound, finally becoming negative. In intervals of the third 
kind, where $” = 4 f'(¢), we discover the same from equation (5), 
so that we may conclude generally that ¢(t) and w(t) are bounded 
for all time in intervals where |4| 2 |y|; by applying the same argu- 
ment to the second equation (4), we establish the conclusion for all 
time without restriction. This of course does not allow us to exclude 
singularities in the solution, but since the right of (4) obviously sat- 
isfies the Lipshitz condition, that is no problem. 

Among the forms of f which this result leaves open to us, we 
shall select the simplest, a cubic polynomial in which for the sake of 
symmetry the square term is omitted; we put 


(4) Sat + KR t*, Rekei Us 


We now discuss the character of the history determined by the equa- 
tions 4. This may be represented by the motion of a particle in a 
plane, whose abscissa and ordinate at a time are respectively the val- 
ues of ¢ and of w which hold at that time; the problem then becomes 
formally a dynamical system with two degrees of freedom to deter- 
mine the path of such a particle. Since the coefficients of ¢’ and y’ in 
equations (4) are always positive, the system is dissipative; since 
it is also bounded, a well-known theorem of dynamics (Birkhoff, 1927, 
Ch. I) allows us to conclude that it will asymptotically approach some 
stable point of equilibrium. The equilibria will be the real solution 
of the algebraic equations obtained by setting ¢’, y’, ¢”, w” equal to 
zero in (4), which are 


Bg? — Berd yp=Ulad + Ks$") 


Boy — Bead? =Ulnp + xy*). (7) 


We shall solve them under the assumption that ¢, = e, = e: the mo- 
tion does not change its character sharply if this is not precisely true, 
as an examination of the perturbation of first order will readily con- 
vince us; and in any case we expect these parameters to be of the 
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same order of magnitude. For numerical applications, a closer ap- 

proximation can always be found. 
The solutions are readily obtained, and, if we exclude those which 

are certainly imaginary, they are in general five in number: 

1. A root ¢, = y. = 0, representing the normal values of the vari- 
ables ; 


2. <A pair of large equal roots 


B(1 —«) | 
Oy eee [62(1 — 2)? —4 py? wa xs}}, 
2 Li k3 We [ks 
Wo1 = her; 
3. A second such pair: 
bCls—a) 1 
Ges == aetas eee ean a [ee a oak 4 w? ki ks ]*, 
ee EL Ks 2u K3 
Woo = do2 5 
4. One pair of large roots, negatives of one another: 
b(1 +e) 1 
al enammmmmmeniie recy ae otek Big (Peo OG 
2 Lt K3 2u K3 
Ws1 = — hai; 
5. A smaller pair of such roots 
b(1 + e) 1 
Serer eae [6?(1 ‘+ €)? — 4 py? ny xs]?, 
3 3 
Ws2 = — dae - 


Some or all of these roots, except of course ¢., y., may happen 
to be imaginary, so that they do not in fact represent equilibria. This 
will be the case for the pairs (2) and (3) if the radical in them is 
negative, a condition which holds if and only if 


Laie Sigs Zieh (8) 


where o = ui ks/f . 

The equilibria corresponding to the roots (4) and (5) will exist 
under an inequality for the radical occurring there which may easily 
be transformed into 


, e>2e-1 
6 hhA Vina et 


For all sufficiently small «, the inequalities (8) and (9) will obvi- 
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ously leave a range for ¢ in which all five equilibria exist. In perhaps 
the most common case, u and x; will be quite small, and f and x, of 
the order of unity, so that o is close to zero, the condition (9) is sat- 
isfied by virture of the sign of ¢; and (8) requires only that ¢ should 
not be in the immediate neighborhood of unity. Still, if uw, x, oF xs 
is quite large or § is small, both inequalities may fail, and all equi- 
libria save the origin cease to exist; this latter is then necessarily 
stable. An organism whose parameters are so related cannot possibly 
exhibit disturbed behavior under any circumstances. 

It is important to investigate the behavior of the solution in the 
neighborhood of the points of equilibrium, first, to determine their 
stability or instability, and second, to throw light on the fluctuations 
of the organism about the normal level which are insufficient to pass 
into permanent disturbance. If 6 = »,, y=», be a point of equilib- 
rium, and if we set 


$=o-n, PHyr-n, 
we shall derive expressions for ¢ and » of the form 


= so 
g(t) =D Ai er 
4=1 
(10) 
= fe 
p(t) =D Ai er? 
1=1 
where the A; are certain linear combinations of the A; , these depend 
on the boundary conditions, and, of the system (4) for this equilib- 
rium, the j; are the roots of the characteristic equation which is 


eae 0 v4 — B & % 
0 aya Br. — 2B &% Va 25) 
1 0 =u 0 
0 1 0 ih 


where 
a, =— [ut+f'()] 
a= — [w+ fi (r2)] 
yy = 2 Bm, — Be: v2 — wh (1) 
yo= pm, ae uf (»), 


provided that all these roots are distinct and that none of them van- 
ishes. This equation may be written in the form 
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(A ae a1) (A =e 2) — Ayala aa 01) ae yo (A ea de) | 


(11) 
+ V1 Ye + Bp? & v2 (v2 —2€%:) a) g 
Since |»,| = |»2| for all the points of equilibrium, we have ao, = a —a., 
and equation (11) is easily soluble, with roots 
a 


MES + fot +2(y. t+ yo) =2V (i + ye)? — Fx)?, = (12) 


in which y = y1 yo + B? € 1 (% — 2 e») and each of the four 4; is ob- 
tained by a different choice for the two doubtful signs in expression 
(12). Now for the equilibria we are considering to be of stable type, 
in the sense that no sufficiently small deviation from that point will 
generate forces which tend to augment it, the condition is that none 
of the 4;, or their real parts, if some of them be complex, shall ex- 
ceed zero. This will be the case here if and only if 


a=—[ut+f'(r)] <0, 
which is always true, 


Yi t+ ye= 386, — 2K — CuK v1 < bern, (13) 


and 
X= (38 ws v1 + Berg + uw — 28 r,) (3 WK; v2 


“| th ky, — 2 ve) = Be 7, OO, — 2 we) On ert 


For the equilibrium at the origin, these reduce to 
eM RG Se LR gee 


so that the origin is always a stable point, as we should expect. For 
the other equilibria, the conditions become extremely cumbersome 
when handled directly; we shall therefore approximate them in the 
following manner. The quantity o? = wv? x, «3/8? may in general be ex- 
pected to be quite small in comparison with unity; we may therefore 
expand in the expressions for the points of equilibrium in powers of 
o°, substitute in (13) and (14), and terminate the expansion at the 
first term that gives a determinate form to the conditions. We dis- 
cover that for the points $a, yoi and ¢51 , Ys1 We May neglect o? itself, 
whereas for the smaller equilibria at go, yo. and $32 , Wee It 1s neces- 
sary to take account of terms in o?; and for one initial value, namely 


€ = 4, we must proceed to terms in ot. The results are as follows: 
the condition (13) is satisfied < 


1) 
2) 
3) 
4) 
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By ¢21 , Yo: unless 3/5 Se <1; 

BY $22, Yoo unless 3 << e [2+ \/10; 
By $a, Ysi if ¢ 2 0, i.e, always; 

By $32, Ws2 if e 20. 


FIGURE 1 
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The inequality (14) is fulfilled 


1) By de, yo unless $< ¢ S1; 

2) By dee, Yoo if e < 4; 

3) By $a, ya if e 20; 

A) By $s, Yo unless V7 —-2<e<1. 


We may obtain a general insight into the qualitative character of 
the motion determined by our equations from the first figure. The de- 
generate algebraic curve consisting of a parabola A together with the 
axis of y divides the plane into the regions for which ¢” is positive and 
those for which it is negative, while the closed curve B together with 
the cubic parabola separate the regions of positive from those of nega- 
ative y”. Clearly, the intersections of these curves will constitute the 
points of equilibrium. In each of these regions, we now have arrows, 
whose slope is the value of y”/¢", which represent roughly the direc- 
tion of motion of a particle placed there with no kinetic energy. The 
presence of kinetic energy and the dissipative forces will modify these 
considerations somewhat, and the initial conditions in (4) will super- 
pose upon these forces an initial velocity with components —f(¢), 
—f(w) toward the origin. 

The difference in character between the stable and the unstable 
equilibria becomes very clear in the diagram. As remarked above, 
these are no periodic orbits in the large, and unless continually dis- 
turbed, the particle will ultimately settle toward one of the equilibria. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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With the view of suggesting a possible mechanism, it is shown that 
sporulation may be synchronous, as with the malarial parasite, given 
the following condition: A metabolite which is essential to the growth 
of the parasite is periodically released from storage organs and dissi- 
pated from the blood stream of the host. 


It is well known that the sporulation of the malarial plasmodium 
exhibits marked synchronism in virtually all species, and that while 
the period varies from species to species, it is in general some mul- 
tiple of 24 hours, with exceptional instances of 12-hour multiples. 
Moreover, numerous experiments (Taliaferro and Taliaferro, 1934; 
Boyd, 1928, 1929 and 1933; Stauber, 1937) have shown that by alter- 
ing the diurnal cycle of the host’s behavior and environment in vari- 
ous ways, the synchronism of the sporulation can be disrupted and 
re-established at a new time in accordance with the new cycle imposed 
upon the host. On the other hand, while it is apparently rather easy 
to disrupt the initial synchronism, the re-establishment is more diffi- 
cult, and there appear to be definite limits beyond which the period 
cannot be changed. 

One can hardly question the conclusion, therefore, that the time 
of sporulation is somehow imposed upon the parasite by the host, but 
that the regulation is not complete, since there are quartan and ter- 
tian malarias as well as quotidian. There seems to be no generally 
accepted theory as to the mechanism involved, although a variety of 
suggestions have been made, based upon attempts to contro] indepen- 
dently this or that feature of the environment. But even if one should 
establish experimentally that, say, the time of sporulation could be 
influenced by the host’s feeding habits, and not by his periods of ac- 
tivity and rest alone, one would still be faced with the task of inter- 
preting this in terms of the physicochemical environment of the para- 
site, to say nothing of the difficulties that would be involved in con- 
clusively establishing such a hypothesis in the first place. 

We wish to consider the problem more directly from the point of 
view of the possible physiological mechanisms that might be capable 
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of exerting such a regulating effect. Many such mechanisms are pos- 
sible a priori, but any one, considered quantitatively, must imply cer- 
tain processes whose presence or absence can be established by direct 
physiological measurements. In this place we present one of these 
possibilities. 

Rashevsky (1940) has shown that a unicellular organism may be 
expected to become unstable upon reaching a certain critical size, on 
account of the diffusion forces resulting from its own metabolic ac- 
tivity. It is indeed possible that, in some cases at least, other forces 
may contribute or may even anticipate the diffusion forces in bring- 
ing about the division, but the diffusion forces are shown to be ade- 
quate and to be effective at the right orders of magnitude. The criti- 
cal radius determined by these forces depends upon the diffusion co- 
efficients, the molecular weights and the metabolic rates for the vari- 
ous metabolites, as well as upon the surface tension and the tempera- 
ture. Also the metabolic rates presumably depend, in part at least, 
upon the supply of the metabolite itself. 

Moreover, the rate of growth of the organism must be dependent 
upon those same parameters, except possibly for surface tension, 
since growth can only occur through the intake of nutrient material 
from the surrounding medium. The simplest picture one can form of 
the entire process is that growth proceeds until the organism reaches 
a critical size, after which division occurs. The synchronism could 
then be due to factors of either or both of two types; factors which 
periodically lower the critical radius, or those which periodically 
yield a marked acceleration in the growth rate. We consider here the 
latter possibility, supposing the critical radius fixed. 

The growth of the cell may be supposed to result from the excess 
of the intake of substance to be consumed over the elimination of 
waste products of the metabolism (Rashevsky, 1940), and the rate 
of intake is dependent, as remarked above, upon the supply. On the 
supposition that the rate of consumption of n metabolites is jointly 
proportional to the supply of these—their concentration within the 
cell—it has been further shown (Householder, 1943) that the various 
consumption rates are all approximately proportional to the available 
concentration of a single one of the metabolites, which thus acts as 
the “limiting factor.” Then, if this limiting metabolite is one whose 
concentration in the blood stream varies sufficiently, the growth rate 
of the parasite, and hence its division time, will be largely controlled 
by this factor. The simplest supposition would be that the metabolite 
is released once during the day from storage organs within the host's 
body, after which it disappears exponentially. The time of release 
might have been established partly by conditioning from the host’s 
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feeding habits, and, though stable for minor deviations in these, yet 
subject to variation as a result of thorough-going disruption. 

Under these circumstances, if we take the concentration of the 
limiting metabolite to be proportional to e-', the relative increase in 
the size of the parasitic cell will be proportional to this, and its rela- 
tive decrease through the process of elimination of wastes will be pro- 
portional to the rate of this elimination, which we suppose constant. 
Rashevsky’s growth equation then has the form 


2rdr/dt = a'e-*t — pre, (1) 
where r represents the radius of the cell. The solution of this is 
r= 1,2 €8! — a(ebt — ect), (2) 
ALLE Th shee FD ORE he dd (3) 
or 
r=—(r7tt+atje*, a=Bp. (4) 


It is enough for our purpose to discuss the case of a< 8. Ifa< f,a 
as defined by (3) is positive, and by a simple change of scale we may 
take a =1. Hence in this case we may write 


72 = (ro? — 1) Bt + et, (5) 


It is natural, furthermore, to take 24 hours as the time unit. Hence 
equation (5) holds for 0 = t = 1 if no division has occurred, after 
which the value 7, must be replaced by a new value, the value of 7 
at t= 1. If this is now taken as a new origin of time, the equation 
holds again in the same form. 

In writing equation (5), 7 was taken as the value of r at the 
time t = 0, i.e., at the time of release of the metabolite postulated. 
As such, 7,” is essentially positive. However, division occurring at 
time ¢t., an equation of the same form will hold for the schizont on 
the time-interval t., 1, provided 7,? is so chosen that the correspond- 
ing r at time ¢, is equal to the radius of the schizont. The 7” so de- 
termined will not necessarily be positive, but if not, this would only 
mean that any parasite, however small at time ¢t = 0, would, by the 
time t = t, , have grown to larger size. 

For discussion it is convenient to introduce, in place of 8, the 


parameter / defined by 
p=a(1-4+.1/A). (6) 
Then in place of equation (5) we may write 
Piste t Lei 1) ee] (7) 
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Also we have | i 3 


2 i 

= aen [ravenna] i (8) 
dt A ied T Bel 4 
Now it is clear that at any time ¢, 7? when regarded as a func- 
tion of 72 is monotonically increasing; moreover, the derivative, and 
hence the growth rate of the parasite, is monotonically decreasing. 
Hence if division does not intervene, two parasites which are initial- 
ly unequal in size will remain so, and in the same sense, but the 

growth rate of the larger one is always less. In fact, unless 


i en © ee Sp (9) 


growth will always be negative, i.e. the size decreases rather than 
increases, And even when relation (9) is satisfied, the right member 
of (8) always vanishes for a single value of t > 0, which may be less 
than unity, according to the value of 72. Hence, in general, except 
for 7,2 small and perhaps even negative, growth ultimately stops, and 
the size of the parasite then decreases until the next release of the 
metabolite. This decrease, however, is slow and need not be detect- 
able. 

If a// is large, say somewhat greater than unity, the second term 
within the brackets in equation (7) becomes negligible for ¢=1. 
This means that the size of the parasite at the moment of the release 
of the metabolite is independent of 7,?, or nearly so, and hence all 
parasites are then of the same size. This independence will be more 
or less complete according to the order of magnitude of a/j, and the 
radius at this time is 


Tr, =ev?, 


Let us suppose then, as we must, that this is below the critical radius. 
Suppose, further, that 


a 22h 
— > log 3 
A a 


Then the curve (7) for 7 =0 rises to a maximum for t < 1, and the 
same is true with any positive r,? satisfying relation (9). If the 
square of the critical radius is less than the maximal r? given by (7) 
when 7,” = e~?, then division will always occur within 24 hours. More- 
over, if 

Tove eae, 


where » is the number of schizonts and ,, the critical radius, then no 
second division can occur. ‘However, if this is not the case, a second 
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division could occur. This would not be in just 12 hours, necessarily, 
but with a deviation of not more than an hour or so, a 12-hour period 
could easily be supposed. In either case the occurrence of synchron- 
ism is obvious. 

If a/d is small, there will be a spread in the radii at the time 
t = 1, but the division may still be synchronous. In fact, the follow- 
ing properties of the curves 


Fa FAG oss 4) (10) 


are the only ones required, and these are easily secured with suitable 
values of the parameters a and /: 


1) f(r? , t) has a maximum in ¢t for t < 1 for at least suffici- 
ently large values of 752; 


2) for each value of t, f(7? , t) is monotonically increasing 
ing; 
3) the condition 
Te — Jlto 5.1) (11) 


is fulfilled for values of 7, not too great. 

Suppose, then, that the critical radius 7, is such that one of the 
curves (10) has 7.2 as its maximum for t= ?t, < 1. Let this be de- 
fined by 7, = 7’, and let 


ry? = f (17,1) > 1 ?. (12) 
Then 
re=f(r'?, tc), 
; 3 13 
0=f,(r'e, te). a, 
Let r”, and 7”, be defined by 
ne yy 2/3 — ii Cee ) i) r) 
14 
vr"? — if (7" 9? y1) . ( ) 
and let 
re yt < 1" 92 "4? > 7.2. (15) 


Then no division can occur for t > ¢., and if division occurs at 
t = t'., then r”,? defines the curve of development of the schizonts, 
and r”, is their radius at t= 1. Conditions (15) insure, respectively, 
that at most one and at least one division occurs in one day, at any 
rate after the first day. 

We can make this clear and establish the synchronism as follows. 
At the start of the second cycle every parasite has a radius 7, satis- 
fying. 7". = ro = 172 80 that every parasite divides during the second 
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cycle and in every one thereafter. Consider, then, any parasite of 
radius 7 at the beginning of the second cycle, yielding schizonts of 
radius 7, at the beginning of the next, 72 at the next, and so on. Be- 
cause of the monotonicity properties of the curves, if 7, < 7, then 
the second division occurs later than the first, hence r, < 7, and the 
third occurs still later, and so sequentially. Likewise, if 7, > 7 , the 
second division occurs earlier than the first, ete. Now if % = 71, 
then necessarily 7, < % , while if 7» = 7"1, 71 > 1%. Hence we can de- 
fine two sequences, one decreasing with initial member 7 = 1,, 
one increasing with initial member 7, = 7", , and no member of the 
second can exceed any member of the first. Hence each sequence has 
a limit. It is easily seen that these limits are identical and define the 
steady-state value of the radius of the parasites at the beginning of 
the cycle. 

We have supposed that each segmenter divides into » exactly 
equal schizonts, and if this were the case, perfect synchronism would 
result in time. Actually, such perfect equality of division is not to be 
expected, and hence neither will perfect synchronism result. Smaller 
schizonts will be more tardy in dividing and some may fail entirely 
during the ensuing cycle. But with our hypothesis, all divisions which 
do occur must do so within a limited time interval, and the distribu- 
tion of frequencies during this interval could be calculated from a 
knowledge of the size-frequencies among the schizonts. Finally, 
whereas we have considered explicitly only a quotidian malaria, it 
should be clear that the same mechanism but with different parame- 
ters »y, a and / could yield equally well a tertian or a quartan period. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of The University of Chicago. 
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A quantitative formulation of Ogle’s induced size effect is carried 
out on the assumption that vertical equality of the retinal images acts 
as a cue for localizing the subjective median plane as opposed to the 
horizontal retinal disparities which serve for localizing the subjective 
frontal plane; that this cue operates independently of such cues as are 
provided by ocular versions, and predominates when the versions are not 
great. The equation which embodies the theoretical predictions involves 
no parameters not directly measurable. 


It is well known, and, from simple consideration of the horopter, 
easily predicted, that a cylindrical size lens placed with axis 90° be- 
fore one eye induces a rotation of the subjective visual space away 
from this eye and toward the other. It is less well known, and not to 
be so predicted, that the lens placed with axis 180° induces a contrary 
rotation which is about equal in amount when the lens magnification 
is not too great. This phenomenon Ogle has designated the “induced 
size effect’ and he has described it in a series of papers listed below. 
It is indicated in these papers that the effect must be somehow due 
to the disparity in the retinal images, but a quantitative formulation 
of this is not made. The purpose of the present paper is to make such 
a formulation and to show how, in at least the “pure” situations, 
Ogle’s experimental results are easily predictable. 

In general, Ogle’s experimental procedure consists in setting up 
a movable plane at a fixed distance from the subject and asking the 
subject to adjust the plane in a position parallel to the frontal plane. 
Ogle notes (1940) that more consistent results are obtained if the 
fusion contours on the plane are “restricted to relatively small areas 
above or below the center of the plane.” In practice these fusion con- 
tours appear to be small circles arranged in two horizontal lines above 
and below the center. 

In the quantitative formulation contemplated, it is assumed that 
an object which forms retinal images that are equal in the vertical 
dimension are judged to lie on the median plane, whereas one which 
forms images that are unequal in the vertical dimension are judged 
to lie on the side of the median toward the eye having the larger 
image. Thus whereas the horizontal disparity of the retinal images 
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provides the primary cue for localizing objects with reference to the 
frontal plane, an entirely different cue, the vertical disparity, is uti- 
lized in localizing the median plane. It is assumed, therefore, that 
Ogle’s subjects were actually locating their subjective median in his 
experiments on the induced size effect, and this seems to be borne out 
by Ogle’s statement concerning the placing of the fusion contours not 
too far from the center of the movable plane which contained them. 
Certainly the two cues which in this case conflict, that provided by the 
horizontal and that provided by the vertical disparities, must both 
operate to some extent. But in placing the fusion contours not too 
far from the center of the plane, the horizontal disparities are mini- 
mized. Also the necessity for placing these above and below the cen- 
ter is plain, since otherwise no vertical disparity would be present. 

To make the necessary calculations, we introduce a coordinate 
system with the x-axis joining the nodal points of the two eyes, the 
positive y-axis extending forward, and the positive z-axis upward in 
the medial plane. It is unnecessary to take account of the small dis- 
tance between the nodal point and the center of rotation of the eye. 
If we introduce as the unit of length half the interocular distance, 
then the left and right eyes have the coordinates (—1,0,0) and 
(+1,0,0). An object whose extremities are at (x,y,0) and (z,y,z) 
subtends at the two eyes angles J, and j, given by 


tan 4, =2[(x $1)? + y?]", (1) 


and if a size lens of magnification wu is placed before the left eye, the 
angle 1’, is given by 
tan 1'5 = pz [(e% +1)? +47]+. (2) 
These angles are equal and the object is judged to be in the median 
plane provided 
17 ec ei ght AW GY Yes Sl (3) 


24 7 2 2 2 
paid -hoy?iass (rte atene\e (4) 
w—- w—1 


Thus as we vary the true frontal plane containing the fixation point, 
the intersection with this of the subjective median describes a circle 


whose center is collinear with the two eyes. The intercepts of this 
circle with the x-axis are at 


(w#1)/(Qu= 1). , 


or provided 


Evidently for 
y > 2u/(w — 1) 
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image equality becomes impossible and the cue fails. 
Experimentally the magnifications used are of only a few per 
cent. If we set 


w=1+u4, (5) 


so that 100u is the percentage magnification, (4) becomes, on neg- 


lecting u? , 
Lee \ ib sean Ne 
U U 


the approximating circles all passing through the origin. For a mag- 
nification of 10%, the error in the radius as calculated this way is 
only 3%. The breakdown mentioned above occurs with the subjective 
median at an angle of 45°, with 


e—y— (1+ 4u)/u. 


If we set 
U 


————— 
Days Se 


then for a given y, 


e 


x 
-=jy[1+eyt+2aytt--], (8) 
y 


and this is the tangent of the angle of deviation of the subjective from 
the true median at the fixation distance y. For small angles we may 
take this as the radian measure of the angle itself. 

We may conclude, therefore, that on the basis of our assump- 
tions, the angle of rotation of the median plane for near fixation is 
approximately proportional to the percentage magnification, to the 
fixation distance, and to the interocular distance, this angle being 
kept less than 45°. At any given fixation distance, however, if the 
magnification is too great, the cue breaks down, since the retinal im- 
ages cannot thus be equalized. We find, in fact, from Ogle’s data that 
the curves graphing angle against magnification, for a given fixation 
distance, actually are nearly linear up to about 6% or so, after which 
they show a decline. The dependence upon the fixation distance is 
borne out also by curves obtained for 20 cm. and for 75 cm. distance. 
It is especially to be noted that the theory contains no undetermined 
parameters but gives absolute measurements. 

The first of the accompanying figures shows the circles (6) 
drawn for magnification up to 10%. The second shows Ogle’s data 
for two subjects (Ogle, 1938), both right and left eyes, compared with 
the theoretical curve. The fixation distance is 40 cm., and this was 
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assumed to be 6 times the interocular distance: y = 12. Actually the 
pupillary distances were given as 62 mm. and 61 mm. in the two cases, 
making y slightly greater than this. These data are not carried to 
the “breakdown” point, and in general this occurs earlier than is pre- 
dicted. In this connection, however, we may note, first, that such an 
angle already substantially exceeds the “small angle” limits for the 
optical formulas employed in the calculation of the magnification, and 
second, that additional cues such as ocular versions, not accounted 
for in the formula, are bound to become increasingly sensible as the 
angle increases. 

This leads up to the general problem of how the various cues 
combine to yield the geometric structure of the perceived visual space. 
Among the purely geometrical cues, disregarding such things as aerial 
perspective, there are at least five: the ocular versions, the vertical 
equality or inequality of the retinal images (on the present theory), 
convergence, the horizontal disparity, and, presumably, accommoda- 
tion. The first two are concerned with the location of the medial 
plane, the second being more sensitive for smaller deviations, but only 
for these. The third and possibly the last serve to localize the gen- 
eral region focussed with reference to the observer. The fourth gives 
depth to this region and gives the cues for small differences. Under 
‘normal” conditions, all of these would cooperate to yield a more or 
less accurate representation of the objective space under examination. 
Under “‘‘abnormal” conditions, whether the abnormality is physiologi- 
cal (ametropia, aniseikonia, strabismus) or is artificially produced 
by instrumentation, the cues may be conflicting and the resulting im- 
pression will be dependent upon the influence of the abnormality up- 
on each of the several cues and upon the relative degrees of dominance 
of the various cues. We are here assuming the complete dominance 
of one or these for small disparities. The present formulation being 
accepted, it now becomes possible to employ simple algebra to calcu- 
late the influence of each of the five cues mentioned when taken alone, 
and thus by artificially varying these among themselves, it is possible 
to determine experimentally the degree to which each contributes to 
the total spatial picture in any given situation. Thus, for example, 
the geometrical and the induced size effects when placed in opposition 
by a spherical size lens appear to combine by simple addition so that 
they just neutralize one another when the magnification is not too 
great. On the other hand, the cues provided by the ocular versions 
do not “take hold” to overcome the conflicting induced effect until a 
minimal critical angle is reached. After this the version cue becomes 


increasingly dominant and overcomes that provided by the retinal 
inequality. 
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The equation for growth by intussusception involves the product of 
the internal concentrations of » metabolites. This product is shown, in 
general, to satisfy a certain algebraic equation of the n-th degree. Ap- 
proximate solutions are exhibited for a somewhat wider class of cases 
than are considered by Rashevsky. 


Rashevsky (1940) derives an equation for the rate of growth of 
a cell by intussusception based upon the assumption that the rate at 
which each metabolite is transformed into the cell substance is jointly 
proportional to the concentrations of all the metabolites within the 
cell. With slow growth, the metabolic rate being regarded as always 
at the steady state, previous derivations give the (average) internal 
concentrations of the various metabolites in terms of their external 
concentrations, the diffusion coefficients and permeabilities, and the 
radius of the cell. Thus there are available the 2n equations 


—M =k We, - (1) 
Cy = Cor + AL» (2) 


for the determination of the internal concentrations ¢; and the con- 
sumption rates —q,. The symbol J é designates the product of all 
the é’s, and the coefficients A; are functions of the parameters men- 
tioned. 

In deriving the growth equation, Rashevsky evades the solution 
of these equations by resorting to an approximation. It is the pur- 
pose of this note to discuss this solution, which requires the solution 
of a single algebraic equation of the n-th degree, and to consider other 
possible approximations. To obtain the single equation we set 


oy, Ax ky ai Co = Col , 


3 
I (€/co) =2, - 
and obtain 
gag=—-h«rlo, (4) 
C:/Cy. = 1 — faze (5) 
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If we multiply together the m equations (5) we have then 
x=1(1—x/a). (6) 


This is the algebraic equation referred to. The ¢ enter the growth 
equation only through the product of all of them, so that it is this 
quantity 2 alone which is needed. 

It is no restriction to suppose that the m substances are so enu- 
merated that 


Then from (5) it is clear that 
eS ay. 
If we graph the line and the curve 
y=2z, y=l(1—«“/a), (7) 


these intersect at the point whose abscissa and ordinate are equal to 
the required value of «. The curve has its y-intercept at 1 and the 
x-intercepts at a; , so that the required intersection lies on the arc of 
the curve between (0,1) and (a,, 0). Consideration of concavity in- 
dicates that the intersection lies between the secant joining these two 
points and the tangent at (0, 1). The intersections of these lines 
with the line (7) are easily found and yield the limits 


1/(1+n/q) 245,1/C1 + 1a), (8) 


where o is the harmonic mean of the a;. These inequalities become 
equalities in the special case when oa, is so small that the other a’s 
do not contribute appreciably to the harmonic mean. From (3) it is 
evident that this case of a small a, occurs when co; is small or A, is 
large in comparison with k;. In this case, that is, the one substance 
alone becomes a “limiting factor’ because of its limited external 
concentration (small ¢,), or its slow diffusion (large A,) in com- 
parison with the k, which determines its consumption rate. 
Rashevsky’s approximation was based essentially upon the as- 
sumption that the a; be, all of them, nearly equal and small. In this 
case, x is small and equal to the common value of the a,. This ap- 
proximation and the one here suggested are exclusive and nearly com- 
plementary. All metabolites with large a’s can evidently be neglected, 
and if the remaining ones have a’s which are small and nearly equal, 
we return to Rashevsky’s case for these which remain. If the small 
a’s are comparable but still fairly different, however, although the 


degree of the equation may be reduced, other methods are required 
for a convenient approximation. 


~~ 
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ON THE ORIGIN OF LIFE 


N. RASHEVSKY 
THE UNIVERSITY OF CHICAGO 


On the basis of the recently proposed new fundamental equation of 
mathematical biophysics, a suggestion is made for a theory of the forma- 
tion of a primitive cell from nonliving material. The discussion includes 
a suggestion for a quantitative formulation of the degree of biological 
organization. It is shown that according to the fundamental equation 
of mathematical biophysics, organization of the nonliving material may 
spontaneously increase under certain conditions, leading to a formation 
of a primitive organism. This process however is a very slow one, re- 
quiring time intervals of several years or even decades. This may ac- 
count for the failure in observing or artificially producing spontaneous 
generation. 


We shall apply here the fundamental equation of mathematical 
biophysics (Rashevsky, 1943c) to the problem of formation of a 
primitive cell from nonliving material. Using the same notations as 
before, we have 


dx; oF oH 
ee (1) 
dt Ox; Chan 


One of the fundamental characteristics of any organism is its 
distinctness from the environment. In the nonliving liquid material 
the natural state is that of a uniform distribution of different solutes 
in a solvent. If a cell or any other organism is to be formed out of 
such a mixture, some of the molecules of this mixture must be seg- 
regated spatially and arranged into definite geometric patterns which 
enable us to perceive that arrangement as a unit distinct from the 
quasi-uniform environment. We may speak of this process as or- 
ganization. 

A great deal of loose talking and writing has been done about 
organization; but, to our knowledge, no one has yet attempted a sim- 
ple quantitative description or definition of it. It is possible that no 
general useful definition can be given at present; in that case, a spe- 
cial definition in this particular instance is indicated, at least as a 
temporary expedient. This is what we propose to do. 

Consider first a purely hypothetical abstract case: namely that 
a cell consists of a unit, composed of certain substances that may in 
general freely spread through the surrounding solution but that ac- 
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tually remain confined to a geometrically clearly limited region of 
space. We may define the degree of organization in this case as fol- 
lows: 

Let 2R, be the linear dimensions of the nonliving medium, for 
simplicity assumed to be roughly spherical in shape. Let the 7) be 
the radius of the cell, which is also assumed to be spherical. Let 


Ry >See. (2) 

When the constituent material of the cell is spread uniformly through- 

out the volume = Re! of the environment, the cell as such is non-ex- 
(3) 


istent. For this case we may set the organization as being zero. If 
on the other hand, all of the constituent material of the cell is con- 
centrated inside a sphere of radius 7), and none outside, we shall 
put the organization of that state as equal to one. Intermediate val- 
ues of organization will be found when some of the concentration of 
the material inside a sphere of radius 7 is greater than outside, but 
when that outside concentration is not zero. 

Denote by x the average concentration of the constituent ma- 
terial in the cell, by y—its average concentration outside. Consider- 
ing the constant amount of the total material in the system, we have, 
because of inequality (2), denoting by C a constant: 


r+ oho =. (3) 


From the above it follows that when « = y, the organization O 
of the system is zero. When y = 0, then O = 1. We may convenient- 
ly define O by: 


O= 2 (4) 


Equation (4) then determines the degree of organization for any 
intermediate stage. 

Thus we illustrate on this simple example how the degree of or- 
ganization can be defined quantitatively. 

In other, more complex cases, different definitions of O may have 
to be given. For instance, we may consider a cell as a sack consisting 
of a membrane composed of a substance m, which holds inside the 
other constituents of the cell, preventing them from spreading out. 
Denoting by x, and ym the concentrations of the substance m in the 
membrane of thickness 6, that is, in the volume 477,26, and in the 
surrounding medium correspondingly, and by x, and y- the average 
concentrations of the other cell constituents, we may define O by: 
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(2% ol Ue)? (aR = Un) 


AS rere ss (5) 
vem Can 
The organization is again zero for a= Yc; Lm = Vow ataisetor 


a uniform mixture. It is equal to one for a completely formed cell. 
The. quadratic, rather than linear, terms are used in order to avoid 
the absurd situation in which a nonuniform distribution of material 
would give O=0. Instead of equation (4) we might also put 


o=(7-") (6) 
a: 


in order to avoid negative value of O. Inasmuch as all the above ex- 
pressions are used only ag illustrations, we shall confine ourselves 
here to expression (4). 

The increase of organization is accompanied by an increase in 
free energy, although the total energy may not change. Changes in 
total energies, connected with changes in organization are in most 
cases rather small. Therefore for simplicity we shall neglect them 
here. But the change in energy transformation may be very large 
when the organization changes. We do not gain much energy by 
crushing a cell, but we either completely stop or strongly reduce many 
energy producing or consuming metabolic reactions. We may there- 
fore consider, at least theoretically, the case in which H does not de- 
pend on O, but F does. Reasons for that dependence may be various. 
Reaction rates may depend in many ways on the concentrations; or 
the confinement of catalists to a small volume may increase their 
efficiency much more rapidly than the increase of their concentration. 
In general we have F = F(O). As an illustration only, we shall con- 
sider here the case 


F=KO, (7) 


K being a constant. 

If F is an increasing function of O, then considering O as the 
coordinate describing the configuration of the system, we see from 
equation (1) that dO/dt will be positive, and that O thus will increase. 

In order to apply equation (1) to our problem we now have to 
compute the viscous force f (dx;/dt). 

This force will be due to the diffusion resistance and may be esti- 
mated in the following manner. 

If the average transport velocity of the molecules is v, then, 
remembering that the average radius of the sphere through which 
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the substance is flowing when converging toward the forming cell 
is $ Ry , we have 
R2yv = : ieee ed 
Pe GO aire 
or 
4R, dy 


3y dt 


(8) 


Considering each molecule as a sphere with an average radius a, 
the force of resistance per molecule is 6z7av, where 7 is the viscosity 
of the solvent. Denoting by M the average molecular weight of the 
molecule, and by N Avogadro’s number, the number of molecules per 


cm? becomes yN/M , and therefore the force of resistance per cm® is 
Sa ES (9) 
M dt 


Hence, because of inequality (2), the total resistance force in the 
system is equal to 


32n?naNR,* dy 


; 10 
3M dt a) 
Putting 
C Ro? 327?naN 
A=—; B=—(>>1); P=—; (11) 
ee Yoh 3M 
we find from equation (3) and (4) that 
A-—(B+1)y 
0 = —___—_. (12) 
A — By 


Solving equation (12) for y, differentiating with respect to time 
and introducing the result into expression (10) we find for the total 
force of viscous resistance the expression 


dO APR, dO 
at (P+ BE RO)’ 


i Introducing expression (7) and (13) into equation (1) we now 
n 


(13) 


PR,*A dO 


(1+B—BO)? dt — (14) 
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APR,* 
B 


and integrating equation (14) with the initial condition O = 0 for 
t = 0, we find 


=a, (15) 


(1+ B)rkt 
ON i (16) 
aBB 
Brkt + 
L+ B 
The time ¢, at which O becomes equal to 1 is given by 
aB 
t, = ——_—_—_-. (17) 
tK(1 + B) 


If, as may be expected, 7 is very small (Rashevsky, 1948c), then 
t, will be very large. Thus the process of formation of a cell may re- 
quire several years or decades. This may account for the failure in 
observing or artificially producing spontaneous generation. 


LITERATURE 


Rashevsky, N. 1948a. “Outline of a New Mathematical Approach to General 
Biology: I.” Bull. Math. Biophysics, 5, 33-47. 

Rashevsky, N. 1943b. “Outline of a New Mathematical Approach to General 
Biology: II.” Bull. Math. Biophysics, 5, 49-64. 

Rashevsky, N. 1943c. “Note on the Hamiltonian Principle in Biology and in 
Physics.” Bull. Math. Biophysics, 5, 65-68. 


* j Fag: i rer ‘: ~*~ ibis Oe ea ate ie, 5; a 


7 ae “ 


- oo ba , , ae : oer 
Ch Ey Aneaa plbeTd pry remus ba thee 


MN ne ee) 


. i } a cere 
lesa) ot actayé Lisdtutmtchtogrs 2 hy See 
Walk 2 ook 


Oiled 


ae 


: Perea Eh a saeco &. inolss 


= 2, 


+ 


INDEX TO VOLUME 5 


INDEX OF AUTHORS 


BLOCH, INGRAM. Some Theoretical Considerations Concerning 
the Interchange of Metabolites Between Capillaries and 
Tissue 


GOMORI, GEORGE, and ALSTON S. HOUSEHOLDER. The Kinetics 
of Enzyme Inactivation 


HOUSEHOLDER, ALTSON S., and GEORGE Gomori. The Kinetics 
of Enzyme Seroeee afl 


HOUSEHOLDER, ALSTON S. Note on Anodal Excitation 


HOUSEHOLDER, ALSTON S. On Synchronous eka with 
Possible Reference to Malarial Parasites 


HOUSEHOLDER, ALSTON §. A Theory of the Induced Size Effect 


HOUSEHOLDER, ALSTON §S. Note on Rashevsky’s Eos for 
Cellular Growth 


LANDAHL, H. D. A Mechanism of Division of a Cell with an Im- 
permeable Membrane : 
LANDAHL, H. D. Studies in the Mathematical Biophysis o of Dis- 
crimination and Conditioning: III .... . : 
LANDAHL, H. D., MCCULLOCH, W. S., and WALTER Pitts. A Sta- 
tistical eaecqaenes of the ocal Calculus of Nervous 
Nets ; She) Pir as 2 learns : 
LETTVIN, JEROME Y., and WALTER Pitts. A Mathematical The- 
ory of the Agecice Psychoses Bam ; : 
MCCULLOCH, WARREN S., and WALTER PITTs. A téateht” Calcu- 
lus of the Ideas Immanent in Nervous Activity ee 
MCCULLOCH, WARREN S., LANDAHL, H. D., and WALTER PITTS. 
A Statistical Consequence of the Logica Calculus of Ner- 
vous Nets F Pee Sng S itOs Sita chs 
PITTS, WALTER. The Linear Theory 0 of Neuron Networks: The 
Dynamic Problem ee re ene : 
PITTs, WALTER, and WARREN S. MCCULLOCH. A Taek: Calcu- 
lus of the Ideas Immanent in Nervous Activity 
PiTTs, WALTER, LANDAHL, H. D., and W. S. MCCULLOCH. A Sta- 
tistical Consequence of the Logical Calculus of Nervous 


Nets io 
171 


PAGE 


149 
155 


161 


75 


103 


135 


139 


115 


135 


23 


115 


135 


172 INDEX OF AUTHORS 
PAGE 
Pitts, WALTER, and JEROME Y. LETTVIN. A Mathematical The- 
ory of the Affective Psychoses 139 
RASHEvskKy, N. Outline of a New Mathematical Approach to 
General Biology: I dest 33 
RASHEVSKY, N. Outline of a New Mathematical Approach to 
General Biology: II aes 49 
RASHEVSKY, N. Note on the Hamiltonian Principle in Biology 
and in Physics . : 65 
RASHEVSKY, N. On the Form of Plants and Animals . .. . 69 
RASHEVSKY, N. On the Theory of Amoeboid Movements . . 95 
RASHEVSKY, N. Mathematical Biophysics of Cell Division . . 99 
RASHEVSKY, N. On the Origin of Life 165 
STANTON, HENRY E. Theory of Cellular Pulsations . . . . 15 


INDEX OF SUBJECTS 


Activity, nervous, 115-133 
Affective psychoses, 189-148 
Amoeboid movements, 95-98 
Animals, form of, 69-73 
Anodal excitation, 91-94 


Biology 
general, 33-47; 49-64 
Hamiltonian principle in, 65-68 


Capillaries, 1-14 
Cell division, 75-81 
Cellular 
pulsations, 15-21 
growth, 161-163 
Conditioning, 103-110 


Discrimination, 103-110 


Effect, induced size, 155-160 
Enzyme inactivation, 83-90 

Equation for cellular growth, 161-163 
Excitation, anodal, 91-94 


Form. 
of plants and animals, 69-73 


Growth 
cellular, 161-163 


Hamiltonian principle in physics and 
biology, 65-68 


Inactivation, kinetics of enzyme, 83-90 

Induced size effect, 155-160 

Interchange of metabolites between cap- 
illaries and tissues, 1-14 


Kinetics of enzyme inactivation, 83-90 


Life, origin of, 165-169 
Logical calculus, 115-133; 135-1387 


Malarial parasites, 149-154 
Membrane, impermeable, 75-81 
Metabolites, interchange of 

between capillaries and tissue, 


Movements, amoeboidal, 95-98 


Nervous 
activity, 115-133 
nets, 135-187 
Neron networks, 28-31 


Origin of life, 165-169 
Parasites, malarial, 149-154 


Plants, form of, 69-73 
Psychoses, affective, 139-148 
Pulsations, cellular, 15-21 


Size, induced, 155-160 

Sporulation, synchronous, 149-154 

Statistical consequences of logical 
apoleks of nervous nets, 135- 


Synchronous sporulation, 149-154 


Tissue, interchange of metabolites 


petwren capillaries and tissue, 


INDEX OF SUBJECTS | 173 


Es BOOK REVIEW 
— . ' PAGE 
_D’ARcY WENTWORTH THOMPSON, On Growth and Form . a 
mee (Reviewed by Gerhardt von Bonin) . £0 3000 3 woe eer Tl 


i 
$ 
ae F 
pars = ‘ 
. Pol 
s - : 
oa a — 
y 
ve < 
fen - 
Ere = 
"v 


is” SREP AAA 
oe at Zakt mpgs i” a) a f 
rin peibget hae , sum mee Sake 
on ‘ F, a ay) Ties t “ as uf ; ; tg E. a 
re vs seg ee Are sien a ile rs sc ek os 
—s e poieii t= | «oN Oe ii Ae as, His: 
are eal bo jautio oy gi ind we hereee 
Bae each) by Ou sear wl ye Shicetion Fe. — ; 
~~ oapeteat oo Tae Sag. ae mis 
i ae me 
ae Po : ber J y 
al b . a= = 
— 


SCOPE OF THE BULLETIN 

1. The Bulletin is devoted to publications of research in Mathe- 
matical Biophysics, as contributing to the physicomathematical foun- 
dations of biology in their most general scope. 

2. Papers published in the Bulletin cover physicomathematical 
theories as well as any other mathematical treatments of biological 
phenomena, with the exception of purely statistical studies. 

3. Mathematical studies in physics or in borderline fields in 
which a direct connection with biological problems is pointed out are 
also accepted. 

4, Emphasis is put upon the mathematical developments, but a 
description and discussion of experimental work falls also within the 
scope of the Bulletin provided that description or discussion is made 
in close connection with mathematical developments contained in the 
same paper. 

5. Outside of the scope of the journal are papers of purely sta- 
tistical nature or papers concerned only with empirical equations. 


PREPARATION OF MANUSCRIPTS 

All manuscripts should be typewritten double space. Equations 
should be numbered on the right consecutively. Do not use such nota- 
tions as “equation 2a” or, “equation 5’” etc. References should be all 
given at the end of the paper, arranged in alphabetic order by authors 
and for each author chronologically, following exactly the style used in 
the Bulletin. In text, reference should be made by giving in paren- 
theses the name of the author followed by the year of publication. In 
case of several publications by the same author in the same year, use 
notations “1940a’’, “1940b’, etc. 

In writing equations, slanted lines should be used wherever pos- 
sible. 

Every paper is to be preceded by a short abstract. 

Drawings should be prepared in a professional manner on white 
paper or tracing cloth following as closely as possible the style used 
in the Bulletin. They should be drawn approximately twice the scale 
of the finished reproduction. Lettering should be made by means of a 
lettering guide. If desired by the authors, drawings may be made at 
the Editorial Office according to author’s instructions, at cost. 
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